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Abstract 

In the instanton models the chiral phase transition is driven by a transition 
from random instanton-antiinstanton hquid and correlated instanton-antiinstanton 
molecules. So far this phenomenon was studied by numerical simulations, while we 
develop alternative semi-analytic approach. For two massless quark flavors, both 
instantons and "molecules" generate specific 4-fermion effective interactions. After 
those are derived, we determine the temperature dependence of the thermodynamic 
quantities, the quark condensate and the fraction of molecules using standard mean 
field method. Using Bethe-Salpeter equation, we calculate T-dependence of mesonic 
correlation functions. 
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1 Introduction 



The nature of the phase transition remains one of the main problems in non- 
perturbative QCD. Not only it is important to understand it from a theoretical 
point of view, but it is also the motivation for current and future heavy ion colli- 
sion experiments. Major efforts are also being made to simulate finite-T QCD on 
the lattice, and many results have clarified to some extent the phase diagram of 
QCD-related theories, see reviews [||, ^. 

At zero temperature the mechanism of spontaneous chiral symmetry breaking 
and even the very existence of most light hadrons is believed to be connected with 
the fermionic (quasi)zero modes generated by instantons. The specific instanton 
liquid model was suggested in |^], and then studied both analytically |Q, |^ and 
numerically 0, ^, ^, 10|. Even the simplest random model (RILM) was found to 
reproduce mesonic and baryonic correlation functions, known from phenomenology 

and the lattice simulations Furthermore, the "instanton liquid" itself, with 
parameters closed to those predicted, were "distilled" from lattice configurations by 
coohng m, H. 

However, it is much less firmly established what happens with instantons at 
finite temperature T. It has been known for a long time that at high T the instanton 
density is suppressed by the Debye-type screening P, H], and it was thought first 
that it may be the reason for chiral symmetry restoration (see e.g. However 



both recent theoretical analysis at low T ||T^ and numerical simulations |1£], show 
that up to Tc the instanton density changes are not significant, and therefore an 
instanton suppression hardly can be an explanation. 

That is why another mechanism was studied in |20|: the phase transition may 
be driven by "pairing" of instantons and antiinstantons into //"molecules" . A two- 
component (or the so called "cocktail") model was proposed, in which at T < Tc 
individual instantons and "molecules" coexist in ensemble, while (for m=0) at T > 
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Tc only molecules survive and the chiral symmetry gets restored. 



In [10 1 the Cocktail Model was studied numerically, by keeping the temperature 
fixed and changing the molecule fraction /. Many mesonic and baryonic correlation 
functions were calculated: they show very dramatic changes as / goes from small 
values to 1. Furthermore, a very interesting hints for survival of some hadrons 
above the phase transition were presented. This scenario of the phase transition 
was been recently confirmed by numerical simulations of the instanton vacuum at 
finite temperature in which both boson and fermion induced interaction between 



the pseudo-particles is taken into account [21|. It was shown that the molecules 
start to form only close to and that they are polarized in color space as well as 
in the Euclidean time direction, as anticipated in |20| . 

In the present paper, we investigate the same two-component (or "cocktail") 
model by different (analytic) methods. The main difference with numerical simu- 
lations mentioned above is that they first calculate quark propagation in a given 
background field (a superposition of instantons), averaging over gauge fields later. 
In our approach we start with the integration over collective variables of instantons, 
deriving effective interaction between quarks. Basically, our next steps are the same 
as in the Bardeen-Cooper-Schreiffer theory of superconductivity or Numbu-Jona- 
Lasinio model of chiral symmetry breaking. Our 4-fermion interaction effective 
interaction is, however, much more complicated, it includes two types of nonlocal 
terms - one from the random instantons and antiinstantons, and the another one 
from the strongly correlated pairs, or molecules. 

The paper is organized as follows: 

In section 2 we calculate the effective interaction ^g// and the overlap matrix 
element T^j, at temperature T, of an instanton and an anti-instanton in the most 
attractive relative orientation whose separation has only a time component (i.e. they 
lie on the Matsubara circle). After that we perform an integration over the 11 quasi- 
zero modes of the molecule to find the corresponding activity (or internal partition 
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function of a molecule) Zmoi- We find out that the equilibrium configuration of 
the // system, in the temperature range around Tc, is when the instanton and the 
anti-instanton are separated by the maximal interval in Matsubara time t = 
they have zero spatial separation and are in the most attractive orientation in color 
group. 

In section 3, we evaluate the partition function, assuming that the fermion 
determinant is a product of the separate instanton, antiinstanton and molecule de- 
terminants. Next we exponentiate them, by performing an Inverse Laplace Transfor- 
mation. Then we integrate over the instanton, antiinstanton and molecule positions 
to obtain an effective four-fermion lagrangian (we study the case of two flavors in 
the chiral limit). 

In section 4, we use this Lagrangian in a Hartree-Fock approximation, lead- 
ing to a theory of free fermions, but with temperature and momentum dependent 
mass. Then we calculate the free energy and minimize it in order to determine our 
parameters such as the molecule fraction /, the quark condensate and the Laplace 
parameters for the liquid and for the molecules. We solve four mean field equations 
and determine the temperature dependence of all parameters. The results clearly 
show the transition from instanton liquid to gas of molecules and the restoration of 
the chiral symmetry. 

In section 5, we calculate mesonic correlation functions, using the Bethe- Salpeter 
equation. We also extract the pion-quark couplings, the pion decay constant /tt and 
its coupling to pseudoscalar current X-j^. All the results naturally respect the Gold- 
stone theorem (for example /tt vanishes at T = Tc) . To compare with the numerical 
simulations for the mesonic correlation functions, wc perform a Fourier transform 
to Euclidean time. Reasonable agreement between them is observed. 
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2 The // Molecules Around 



This section is relatively independent of the rest of the paper in the sense that it deals 
with properties of the single II molecule. It consists of two steps: the first is a new 
calculation of the // interaction in the relevant configuration (so to say, the quantum 
mechanics of it); while the second deals with its internal statistical mechanics. In 
the next sections we will proceed to the statistical ensemble of multiple instantons 
around the critical temperature. 

The // interaction and the overlap matrix elements of the fermionic zero modes, 



both at zero and at finite temperature, were studied in detail in |22, 23|. Unfortu- 
nately these studies were not complete in the sense that they did not include our 
region of interest - the configurations in which the pseudo-particles lie around the 
opposite ends of a diameter on the Matsubara circle. 

The best definition of the classical II interaction Sint = {S — 2Sq), where Sq 
is the single instanton action, is given by the so-called streamline configurations 
p4| . Unfortunately, at finite temperature, the conformal symmetry that allows to 
find the streamline is missing. Therefore, one has to use some ansatz for the // 
configuration which has a natural extension to non-zero T. We are going to use the 
ratio ansatz ^ 

^'^^''>-—g n, + nj-i ' 

where 11/ = n(x — Zf), Hj = Il{x — zj) and 

1 '^T^p'^ sinh(27rrT) 
^ r cosh(27rrT) -cos(27rtT)' ^' 

which is free from the artifacts of the sum ansatz and provides a reasonable repulsive 
core. The classical bosonic interaction is just obtained from this expression, by 
numerical calculation of the classical action. 

The quark-induced interaction is described in general by fermionic zero-mode 
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overlap matrix element: 

Tjj = J dt(fx(j)\{x - zi) ^(piix - zj), (3) 



In this case we use the simplest sum ansatz [ 25 1 , because in this case one can use an 
ordinary derivative instead of a covariant one. As it was shown in [^], the results 
in this case differ insignificantly from the ones obtained by the ratio ansatz. 

The 0(3) symmetry of the // configuration along the Matsubara circle allows 
us to reduce the above integral and the integral for the bosonic part of the action 
to a two dimensional one, and we perform the integration numerically. Our results 



coincide with the ones from [ 22 , in the region where their formulas are valid (for 
Matsubara time separations that are not close to L/2 = The results for Tjj 
are shown in Fig. 1 (a), and the combined e~^'-"''Tjj in Fig. 1 (b) (the dotted line). 

The second step is to include the "internal entropy" , or multiplicity of different 
// configurations into account. One has to take the integral 

ZmoiiT) = J d^'^l T^f ^e-^'- (4) 

over all collective variables. These are the time separation r, the three dimensional 
spatial separation and, for SU{3), seven relative orientation angles|^. Most of the 
11-dimensional integral, however, can be done by a saddle-point approximation. 
Both Sint and Tjj decay exponentially in the spatial separation directions. The 
second derivatives of Sint over orientation angles are big, because they are multiplied 
by Svr^/gi^ = 5o ^ 10 >> 1. In contrast to that, the angular dependence of 
— 2Nflii(Tjj) is rather weak, so in angular integrals we treat Tjj as a constant. 
Finally, we find that the saddle point approximation is not good at all for the 
integral along the time separation r 0. That is why we perform this last integration 
explicitly. 

^ Note that rotations in tlie direction of Ag does not ciiange the // configuration, so that we 

have only 7 orientation parameters. The measure of integration is described in [pot- 

E.g. the second derivative of Sint — 2Nf ln{Tjj along this direction for a time separation 



In Fig. 1 (b), wc show the time dependence of the integrand Z'^^iij) = 
J T^jC^"^'"', (the integration over ah variables but the Matsubara time), at Tg, 
compared to the exponent of the classical and quark- induced interaction exp —{Sint— 
2Nf ln(rf j)}, which has a maximum at about t = L/A and a minimum at r = L/2. 
One can see that when the multiplicity of phase space is included, the II molecules 
are distributed around r = 1/2T. The corresponding time dependence of Tjj at Tc, 
which is is shown in Fig. 1(a), is rather weak around the maximum of Z'. This 
allows us to use its central value for all molecules that participate in our "cocktail 
model" . After the last integral is done, one gets the absolute value of the statistical 
sum for molecules Z^oi. The temperature dependence of Z^ol and Tjj is shown in 
Figs. 1 (c), (d). 

How accurate is our saddle point integration? In the next section we will find 
that in order to set the phase transition at Tc = 150MeV one needs Zj^oi about 1.8 
larger: it means only a 5% — 10% error per integration over each of the // collective 
parameters. However, we think real uncertainty in Z^oi is larger, and it is about 
one order of magnitude (as one can judge from including — 4 In (Tjj) in the second 
derivative over the orientation parameters, etc). Nevertheless we think its tempera- 
ture dependence is evaluated reliably, and we use it below in our calculation of the 
thermodynamical properties of the instanton ensemble , the correlation functions 
and the pion coupling constants. 

3 Effective 4-Fermion Interaction 

Let us start with a list of assumptions and simplifications, defining the model and 

the particular parametrization used. We emphasize, that most of them are possible 

only because we are primarily interested in a comparatively narrow interval of tem- 
of L/2, is negative for temperatures T < 230MeV, although exphcit calculations show that this 
configuration is a maximum for z and not a minimum. 
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perature T around Tc, in which significant structural changes (the phase transition) 
takes place. 

• We fix the total pseudo-particle density to be a constant for T up to Tc = 
150MeV, with only composition / changing. We take nt = Ifm^^ specifically, 
lacking more accurate number. (None of the results change qualitatively if, 
say, it is modified by a factor of 2.) 

• Lattice data in |jl^ indicate, that the instanton radius doesn't change much 
with T till well above the phase transition. So we also fix it at p = 0.34/m. 
(However we are going to express all dimensional quantities in units of p and 
its powers anyway.) 

• In the previous section we have shown that the dominant configuration for 
the // molecules is the one, in which the distance between their centers in the 
Euclidean time direction is equal to "half-box" R = L/2 = 1/2T, the spatial 
distance is 0, and they have the most attractive relative orientation in color 
space. Although the molecule distribution has some span around L/2 in the 
time direction, T^j changes little, which allows us to use one value of it for all 
molecules. 

• Tc and 5*0 are related by the requirement that all unpaired instantons dis- 
appear at Tc- This condition can be satisfied only in a window around 
Tc = Ih^Mev. We use the following values: Tc = ISOMeF, = 9.5. How- 
ever there is a big uncertainty in this relation, due to the uncertainty of the 
calculated value of the molecule activity z. 

• We are going to ignore all bosonic and fermionic interactions between random 
instantons and among // molecules. 

Now we proceed to calculation of the effective quark interaction. For one in- 
stanton it is the well known 't Hooft effective interaction. It can be derived from 
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the zero-mode part of the quark propagator : 

at \ c f \ , (t>i{x - zi)(j)\{y - zi) 

S[x,y) = SNZM{x,y) ^ — , (5) 

—im 

where (j){x) is the quark zero mode in the field of an instanton with center at 
zj and SNZA4{x,y) represent contributions of non-zero modes, (following Q we 
shall approximate these contributions by the free quark propagator So). For Nf 
light fermions, one should take the corresponding power of S{x,y), multiply it with 
the instanton probability (P~ m^i) and then cancel m in both, resulting in the 
interaction which is finite in the chiral limit (m=0). 



Similar procedure can be repeated for an isolated // molecule |21]. In the chiral 
limit quarks have the following propagator: 

^ ^ ^ ^ (t}i{m){x - {zm- L/A))(t)\ Ay - {Zm + L/A)) 
Sm[x,y) = So[x,y) H — h 

<Plim){x - (Zm + L/A))(l)\^^^{y - {Zm - L/A)) 

^11 

where T^j was defined in (2). 

For Nf = 2, the following fermionic path integral generates the above propaga- 
tors: 

+ ■ ■ 'm- J— I m=l 

where 

2 
2 



e. = J dzjd^j n ( / d^x^]{x)i ^(t>j{x - ZJ) J d^^iy - zj)i ^^f{y)), (9) 

dZmd^m [ n (^// + / d*xi)Yx)i ^(l)l{m) {x - Zra + -)x 
f=l 

/LP L 
'^^y^f(m)^^ 4)^ ?>^f{v) + J d^xi;^f{x)i ^4>^^){x -Zm-^)x 



L f L 

d'^y^\(m)(y -^m + -T)i ^1pf{y) / d^x'ipj{x')i ^(/)J(„) (x' - - -)x 



Note, that there are no odd terms in Tjj, because of the integration over 17^- Here 
Zi and rife are the collective coordinates of an // molecule, c = Zmoi/Tjj- and Cp is 
the single instanton partition function integrated over the instanton radius p. 
The convergence of the latter is due to the non-perturbative g{p), which tends to a 



constant for large p |14 



Now, applying inverse Laplace Transformation, we get the following partition 
function: 



Z = const J dp+dp-dPmexp{-{N+ + 1) log(/3+/cp) - (7V_ + 1) log(/?_/cp) 

-(Nm + 1) log{Pm/{clc)) + / d^x{lpU - 0+9+ + 0^6^ + PM}. (H) 



To evaluate the effective interaction 9i terms we go to momentum space and do the 
integrations over the orientation angles and the center coordinates of the instantons, 
antiinstantons and molecule^. The Fourier transforms of the fermion zero modes 
at finite temperature are j2g[|: 



a(k,u;„) = (k2+a;2)(A2 + i?2), (12) 



1 f°° , 9 , /"^^"^ , f cosikr) smikr) \ , ,„i/9 
A = / A-rrr^dr dt [ — — ^ -j-^ ) cos{ujnt)Il^' x 



2ttp Jo Jo \ kr k'^r'^ 

^ / (n-1) cosTTtr ^^^^ 

n coshvrrT/ ' 



1 1"°° . 2 , f^/^ , sm(kr) , , ,/2^ / n-1) cosvrtr \ , 
27rpio io kr ^ ' V H coshvrrry ' 



^ Let us again remind that for the molecules we assume complete polarization both in coordinate 
and color space. 
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where 11 is the potential (2) that also appears in the finite temperature zero modes 
(j)^ in coordinate space [30|]: 



^nv^|,fffl^^).«, (15) 



2Trp \ n coshvrrT^ 

and is the standard constant spinor couphng spinor and Dirac indexes. Com- 
bining ah fermion terms together, one obtains the fohowing nonlocal four-fermion 
interaction 0: 



4/ 



n 



1 



i27r)^5\ki -k2 + k3- k4)^a{ki)a{k2)x 



0(^3)0(^:4) I (/?+ + P-)^^ (V'V^ V)^ - (V'^T-fl + octet terms 

^ ((^VaV)^ + (V'^a^sV')^) - ((?/'"^7-a7oV')^ + (V'^'^aTOTsV') 



4Af2 



+ 



^(^^ToTsV')^ + octet terms |, (16) 



where = {i, r), while = (1, r). The four fermion interaction is given in a Fierz 
symmetric form, which means that all permutations of these terms are included. We 
have shown explicitly only the color singlet terms. The last square brackets is the 



"molecular interaction" derived in [21|. There are no two- fermion terms, because 
the integral over the orientation of the molecule in color space is zero. 

One can view this Lagrangian is a variant of the (nonlocal) Nambu- Jona-Lasinio 
model at finite temperature, with particular coefficients and with a natural cutoff 
given by the nonlocality of the vertices, at a scale 1/p. We can now treat this 



interaction analogously to refs. [31, |3^ for the NJL model. 



^AU expressions throughout the text should be understood as given at finite temperature even 
if we didn't write them exphcitly in this form for simpficity of the notations, e.g. J d*k/{2Tv)^ 
meansTX;r=-cx,/ d^^/i^^f' etc. 
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4 The Mean Field Approximation and Thermodynam- 
ics 



To calculate the thermodynamic properties of the system, we first bosonize the 
fermionic action, making Hubbard-Stratonovitch transformation, 

^n«(^0(^V)' ^ 2^a{ki)a{k2)i}\ki)ij{k2)VA{27:)H\k^ - ki)Q- 

Vi{27rf5^{ki - k2)5\k3 - kA)Q^. (17) 

Then we consider the bosonic field Q as a constant, to be determined from free 
energy minimization. As there is no net topological charge, /3+ = /3_ = /?,. We also 
define / = ^^,N^ = = ^-Y-ntotVA, = ^ntotVi- After integration over the 
fermion degrees of freedom we have the following free energy^: 

^ = -TfjPm + (1 - fhtot log(^) + Ltot log(|^) + ^iP + Wm)- 

4iVc/ -0^,\og{k^ + M\k)\ (18) 

where one can identify the chemical potentials for the "liquid" component log(^) 
and for the "molecular" one log(^c),. The last term corresponding to a ga^ of 
massive quarks with the (instanton-induced) momentum-dependent mass 

M(k) = °(*)W + ^/^^> (IS) 



^ Note that, our free energy is not complete in the following sense: the gauge fields are assumed 

to be only a superposition of instantons, while all excitations of the gluonic degrees of freedom are 

excluded. In principle, there should be a gluonic term (similar to the the last quark term in the 

action above) , which eventually (at high T) will lead to the perturbative gluonic part of the thermal 

energy. We do not include it because (i) Glueballs are much heavier than mesons and constituent 

quarks, and are not excited at T ~ Tc = 150AfeV. (ii) We mainly use the free energy in order 

to determine parameters such as / and the quark- related quantities Q, The missing term 

describing gluonic excitations can hardly depend on them. 

®We remind that the instanton vacuum has no confinement, so quarks just change their mass 

in the transition. 
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We shall determine this mass from a self-consistency condition, represented graphi- 
cally in Fig. 2 (a), where the vertex K is a sum of the instanton and the molecular 
interactions (Fig. 2 (b)). 

Our next step is minimization of AF with respect to /, (3m, P, Q- There are four 
equations: 

(20) 

rj.2 f^tot^^^ f d'k 2M{k) Qa{k) 
(1 - f)ntot 

J (27r)4 fe2 + M2'(^) 4]^! ~ STVf " ^^^^ 

One can notice, that the first equation is the condition for chemical equilib- 
rium between the random and molecular components. These four equations can be 
reduced to one "gap" equation, that has to be solved numerically: 

1 f d^k {l + 2cp)a^{k) 



(3m 




= 0, 


d^k 2M{k) Qa{k) 




= 0, 


(27r)4/c2 + M2(A;) 4N^ 


4iV2 


d'^k 2M{k) Qa{k) 


Q' 


= 0, 


{2itY k'^ + M\k) AN^ 


87V2 


f d^k 2M(k) a{k) 


2Q 


= 0. 


1 (27r)4 A;2 + M2(A;) 47V2 


87V2 



+4iVc , :'r..w -ik?^. = 0, (22) 



N, J (27r)4 fc2 + (1 + 2c/3)2(l - f )ntotaHk)/{2N^) 



(24) 



with: 



1 



/?=[-(!- f)ntot + - f?<t + f/i^c)ntotTfj]-^. (25) 

At Tc, all unpaired instantons disappear (/ — 1) and the chiral symmetry is 
restored. The condition / = 1 is satisfied on a line in the T, Sq plane . We have 
chosen = 150MeV = 0.26/p and = 9.5. Of course, when / = 1 or / = 0, 

the saddle point approximation for the /? integrals is no longer valid, so our model 
is restricted in the region < / < 1, where the saddle point method is justified in 
the thermodynamic limit V4 —>■ 00. We see that when f 1, (3m = ^/2ntot/{iTjj)'^ , 



(3 = ntot/i2c{iTjj)'^) does not tend to 0, so that the C/a(1) symmetry remains 
broken when T Tc. This can be seen also in our discussion of the correlation 
functions in section 4. 
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Our results are presented in Fig. 3, where the molecule fraction /, the effective 
quark mass M(0, vrT), the quark condensate {"ipil^) and the energy density e are 
presented. Note, that the local (i/^ip) has one power of the formfactor a less under the 
momentum integral, compared to the nonlocal mean field Q. {ipip) is normalized to 
the phenomenological value of the quark condensate at T = 0. We see that in about 
20MeV the fraction of the molecules drops to about half, the quark condensate 
almost reaches its phenomenological value at T = 0, but the constituent quark 
mass, which is defined as the value of the momentum dependant mass from the 
gap equation at zero spatial momentum and Matsubara frequency vrT, is still low - 
about 1/3 of its phenomenological value at T = 0. (The reader should be warned 
that one can hardly compare the mass value obtained in the present calculation 
with the usual constituent quark mass M{k = 0, u; = 0) ~ SbOMeV because the 
minimal energy possible u ~ vrTc is rather large.) In addition, the finite temperature 
formfactors have a maximum at spatial momentum \<.max 7^ 0, so we show with a 
dotted line the temperature dependence of the mass defined at k = {'kmax,'^T). 

This behaviour of the thermodynamic quantities looks similar to the one ob- 



tained in |2£] or in the NJL models [^], and both are governed by similar gap 
equations for the constituent quark mass or the quark condensate. However the 



physics of the transition in our case is quite different. In |32| the transition is gov- 
erned by the suppression of the quark condensate by thermal excitations of quarks. 
In ||2^ there is an additional thermal suppression of the instanton density, which 
results in additional suppression of the effective coupling constant. In our case, the 
main effect is the reorganization of the instanton vacuum: it leads to reorganization 
of the Lagrangian by itself. This can be seen by the fact that the mean field Q 
(a nonlocal analogue of the quark condensate), is proportional to square root of 
the density of the random component 1 — /. Our mechanism makes the transition 
region narrower and the transition temperature Tc lower. 

In Fig. 4 (a) we show the pressure p = —AF (the solid line), together with its 
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different components: pressure of free massless fermions ANc J log(A;^) (open 
triangles), the deviation from it, due to the effective mass ANc J ^og{^--^-^-^) 
(open squares), the contribution from the condensate —-^^{13 + 2j3m) (stars), the 
contributions from the molecules T'jj(3m — ■^'i-toi log(|z^) (black squares), and the 
instanton liquid —(1 — f)ntot log(^) (black triangles). For comparison, we also show 
that the pressure of a pionic gas p = 0.987T^ (the dotted line) is significantly smaller 
than any of the components under consideration, and thus unimportant. 

In Fig. 4.(b) the energy density e below the phase transition point is shown 
to be actually directly proportional to the molecule fraction /. This correlates 
well with the observation made in ||2l|, that (unlike the individual instantons) the 
molecules have a net positive energy, even in the classical approximation. 



5 Mesonic Correlation Functions 

Our last step is investigation of the effect of the 4-fermion effective interaction 
on mesonic spectra at T ~ T^. Using Bethe-Salpeter equation one may calculate 



mesonic corelation functions, similarly to what was done in |31, Our major 

advantage is that we naturally have a non-local vertices, which provide an ultraviolet 
cutoff. We start with the two-body interaction kernel K, which is given by the four- 
fermion terms in the effective lagrangian. Then we have the BS equation for the 
quark-antiquark T matrix(Fig. 2 c) (See also Fig. 2 (c)): 

T{q)=lC + iTr J ^^a{p+^)a{p-^)[JCSF{p+^)T{q)SF{p-^)]. (26) 

The trace is taken over the Dirac, flavour and color matrices. For the colorless meson 
channels we need only the color singlet terms in the lagrangian. Then using the 
symmetries of the Matsubara space-time, we can decompose T and K into covariant 
structures. 
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(g4g3|/C|g2gi) = E^a 



34 



(27) 
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where Fa denotes Dirac tensors. 

Because at T 7^ there is less symmetry than at zero temperature, we have more 
structures. However our particular action (15) contains only scalar and pseudoscalar 
terms and the time oriented vector and axial terms, produced by 70 and7o75, there- 
fore only the following coefficients are non-zero: 



Kh = 



^[2/3^ + 2(<5io-^)/3], 



V 



4iV2 
1 

2N} 



1 



2iV2 



(1 - 45io)/?r, 



(28) 



If we define the loop integral: 



d^p 



JWi) = J (i^y4«(p + i)"(p - i)[ra ySf(p + |)r/3ySF(p - 1)], (29) 

then the solution of the BS equation (in matrix notation) is: 

T =[1- JT]-^K. (30) 

To get the mesonic correlation function for the corresponding channel, we have to 
multiply T from the left and the right with two loop integrals J (those have only 
two formfactors ^/a, instead of four, because the correlation function is defined 
for point-like currents). For all channels, we have separated equations except for 
the pseudoscalar and axial vector channels which mix both for isospin 1 (vr) and 
(rj). In Fig. 5 we show the Fourier transforms of the correlation functions for 
Euclidean time separation of the currents, normalized to the free correlation function 
at finite temperature. We show them in 6 steps with AT 9MeV, starting from 
T = 150MeV and ending with T ^ 105MeV. 
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The most striking feature is the strong attraction in the pion channel, which 
remains robust and does not disappear at Tc- This is due mainly to the residual 't 
Hooft interaction (induced by the propagating quarks) and in a lesser extent to the 
attractive character of the molecular interaction in this channel (triangles). This 
feature is in agreement with the lattice simulations around , which also show 

a strong pion signal beyond Tc- Futhermore, at Tc our pion signal is similar to the 
one from the numerical simulations of the interacting instanton liquid [^] , although 
stronger]]. 

The results are in agreement with all chiral theorems, and they clearly show 
restoration of the chiral symmetry at Tc- In particular, at T = Tc the vr pseudoscalar 
correlation function {PP^) coincides with the a scalar correlation function {SS^). 
Another feature in agreement with the chiral symmetry is that the pion is decoupled 
from the axial current (/^r ^ 0) at T — > Tc, but remains coupled to the pseudoscalar 
one (Att 7^ 0) even in the pure molecular vacuum (f=l). 

An open theoretical issue debated in the literature is how strongly the f/yi(l) 



chiral symmetry is violated at T ~ Tc, see e.g. 36, 37 1. We remind that 

this symmetry is strongly broken by random liquid at T = 0, but it is respected 
by the new term in the Lagrangian due to "molecules" which we derived above. 
Although molecules are prevailing above Tc the 't Hooft interaction does not dis- 
appear completely at any T. One way to explain it, is to say that external quarks 
induce additional instantons, absent in vacuum. 

The way to measure [/yi(l)-violating effects is to calculate the correlation func- 
tion for r] (the isoscalar pseudoscalar) or the isovector scalar, to be referred to as 6, 
and compare their properties with the pion or sigma ones. These two channels are 
[/yl(l) partners of a and vr, and if this symmetry gets approximately restored, those 
should converge. 



^ The differences rnigfit be due to tfie fact that we are considering the chiral hmit, while in 
the quarks have a nonzero mass, which leads to the smearing or all signals. 
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As shown in in the random instanton vacuum the correlation functions in 
the rj, 6 channels, display so strong repulsive interaction that they become negative^. 
Our results show that although the rj, 5 correlators increase when T approaches the 
phase transition point, the U{1) restoration still does not happen in our model. For 
comparison, we have plotted in Fig. 5 (triangles) the same correlation functions at 
Tc if only the U (1) symmetric molecule-induced interaction is included. 

Because we obtain the correlation functions in for Euclidean momenta by nu- 
merical procedure, we cannot analytically continue them and find the pole meson 
masses. However, the tt mass is guaranteed by the Goldstone theorem to be 0. We 
can therefore derive the pion decay constants /^r and Xt^. The quark T matrix in 
the pseudoscalar, 1 = 1 channel has the following form near the pion pole (see Fig. 
2 (d)): 

rr I _ [9p{il5T3) + ^^(-^ToTsTa)] O bp(n5T3) + (^70757-3)] . 

-'7r|qo=0,q^O - o , ^2 ^■^^) 



Comparing with (30), we get: 



gp 



lKl{l-KlJf,^) 



9A = gp- 



dD^/dq^ 
7^3 t33 \ 



go=0,q-»0 



go=0,q^O 



(32) 

(33) 



'1 i^^J^o^o) 

where = det(l - K^J^) = (1 - KlJfp){l - K\jfj^) - KlK\{jfj^f. Using the 
definitions: 

K^ij = (0|V;(0)i75yV(0)|7r,(g)), (35) 

and calculating a simple loop diagram (Fig 2. (d)), we get: 

= fppgp, (36) 
®This is an artefact of the "randomness" of the instanton hquid related to too strong fluctuations 
of the topological charge: it disappears in the interacting instanton liquid. 



18 



fn = ifl>Ao/q'gp + floAo9A/q')Uo=0,ci-^0, (37) 
fl = iTpA^/Q'9p)\qo=o,ci^o, (38) 

where the tilde indicates that the loop integral is with two factors of y/a only, and 
because of the lack of 0(4) symmetry, there are two /,r's, coupled to the time and 
the spatial components of the axial current. The results are shown in Fig. 6. and 
the two /tt's go to at Tc as required by the restoration of the chiral symmetry. Ayr 
however, remains finite. This means that the pion (and also his chiral partner, a), 
survive the phase transition! This conclusion is consistent with numerical evidences 
obtained from the calculation of the Euclidean correlation functions in the time 



direction [10|, but there Ayr although remaining finite, decreases towards Tc, while 



in our calculation Att slightly increases. This difference, as we have mentioned 



before, might be due to the non-zero quark current masses in |10|. 

Furthermore, in our approach we have found a signal of attractive interaction 
also in the vector (p) channel. It is seen as an additional maximum in the correlators, 
shown in Fig. 5 (i), which is the same as the one, generated by the molecule-induced 
interaction (triangles). This signal was not observed in [^0. 



6 Summary 

In this paper we have studied individual instanton-antiinstanton molecules, in much 
greater details than it was done before. We have determined classical bosonic inter- 
action in the ratio ansatz and quark-induced interaction. Then we have performed 
'One possible explanation is related to the different treatment of the non-zero mode part of 
the quark propagator. In an additional repulsive interaction is included. The same authors 
have done simulations with the free propagator, that we use, for the non-zero mode of the quark 
propagator and they have noticed an attractive correlation function in the p channel. However some 
other contributions are not included in both approaches (e.g. the confinement, which provides an 
attractive interaction). So the question, whether p "melts" below or at Tc, remains open. 
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11-dimensional integration over collective variables, 10 in the saddle-point approxi- 
mation and the last one explicitly. 

The results obtained, have been used in studies of molecule formation in the 
ensemble, at temperatures close to chiral restoration point T ~ Tc. We have used a 
two-component (or "cocktail" ) model, with contributions from uncorrelated instan- 
ton liquid and polarized // molecules, and have confirmed that chiral restoration is 
driven by formation of the instanton-antiinstanton molecules. 

Both random and molecular components have generated an effective 4-fermion 
interaction. With standard mean-field methods we have derived semi-analytically 
the thermodynamics of the system. The basic conclusion is that there is a rapid 
temperature dependence of the fraction of molecules f, see Fig.3(a), which jumps 
from / ~ 0.5 at T = 0.7Tc to 1 at Tc- The corresponding jump in energy density is 
strongly correlated with it (see Fig.4(b)), confirming the idea suggested in |2l| that 
formation of "molecules" is the major reason of rapid growth of the energy density 
around the phase transition point. 

We have also calculated the mesonic correlation functions, using Bethe-Salpeter 
type equation. We have found that the pion is decoupled from the axial current at 
Tc, as it should, but remains coupled to pseudoscalar one even in pure molecular 



vacuum (f=l). A strong signal, similar to that reported numerically in |10|, has 
been found, indicating that pion may survive the phase transition as a bound state. 
Sigma meson follows the pattern of chiral symmetry restoration, joining the pion. 
The "repulsive" channels r/, 5 show increase of the correlators, but they do not 
exactly join the pion and the sigma ones, so U{1)a symmetry remains broken. 
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figure captions 



figure 1 (a) The time dependence of Z'^^^ = J d^^QTjje~^^"* , (the integration over 
the Matsubara time is missed) at Tc, (b) the time dependence of Tjj, at Tc, (c) the 
temperature dependence of Zj^oi, (d) the temperature dependence of Tjj. 

figure 2 (a) The Hartree-Fock (self-consistency) equation for the quark propagator, 
(b) the four fermion interaction is a sum of the 't Hooft vertex and the molecular 
vertex, (c) The Bethe-Salpeter equation for the T matrix in the mesonic channels, 
(d) The pion T matrix near the massless pion pole, (e) The diagram for the pion 
coupling to the axial current. 

figure 3 (a) The molecule fraction /, (b) the effective quark mass M(0,7rr) (solid 
line) and M(k„jaa,, vrT) (dotted line), (c) the energy density e and (d) the quark 
condensate {tpip), normalized to the phenomenological value at zero temperature as 
a function of the temperature. All dimensional quantities are in units of the inverse 
instanton radius 1/p = bSOMeV. 

figure 4 (a) The pressure p = —AF, and the different components of it, as a function 
of the temperature. 

The pressure is shown by a solid line, and the different contributions to it are as 
follows: pressure of free massless fermions iNcJ -^^loglk"^) (open triangles), the 
deviation from it, due to the effective mass AN^ J log(^-^p-^) (open squares), 
the contribution from the condensate — ^^(/3 + 2/3^) (stars), the contributions 
from the molecules T'^jPm — 2'^totlog(^) (black squares), and the instanton liquid 
— (1 — f)ntot log(^) (black triangles). For comparison, the pressure of a pion gas is 
shown by a dotted line. 
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The the dependence of the energy density e on the molecule fraction / (fig. 4b) . 



figure 5 Mesonic correlation functions for Euclidean time separation of the cur- 
rents, normalized to the free correlation function at finite temperature. At r = 
and r = l/T they are 1, but we have subtracted 1 for convenience, r is in units of 
p. These are the isospin 1 pseudoscalar-pseudoscalar (a), pseudoscalar-axial (b) and 
axial-axial (c) correlators, the isospin pseudoscalar-pseudoscalar (d), pseudoscalar- 
axial (e) and axial-axial (f) correlators, the scalar-scalar correlators for isospin 1 (g) 
and (h), and the vector- vector correlators,which is the same for isospin 1 and (i). 
There are 6 different temperatures, starting from T k. 150MeF and ending with 
T 105Mey. For comparison, the correlation functions of a purely "molecular" 
vacuum at Tc are plotted with triangles. 

figure 6 The effective pion-quark couplings gp and gA (a) and the pion constants 
, and /* (b) as a function of the temperature. All dimensional quantities are 
in units of the inverse instanton radius 1/p = 580MeF. 
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